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AN INVESTIGATION OF THE FEASIBILITY OF 

OBTAINING LATERAL STABILITY DERIVATIVES 

FOR A LINEAR AIRCRAFT BY MATCHING ANALOG 

COMPUTER TRANSIENT RESPONSES TO FLIGHT 
TEST DATA 


SUMMARY 


Using a North American NAVION as a flight test vehicle, and an 
electronic analog computer as a simulator, a study was made of the 
feasibility of obtaining certain lateral stability derivatives by matching 
the transient responses of an analog computer on which the lateral 
operational equations of motion for a NAVION aircraft were set up, to 
the actual airplane lateral response in sideslip angle, roll rate, and 
yaw rate. 

The matching process was achieved by supplying time histories 
of the actual recorded control motions as forcing function inputs to the 
computer, and systematically varying those coefficients of the analog 
which corresponded to the stability derivatives of the dynamic equations 
of motion. 

The method was demonstrated to be a practical procedure for 
obtaining the lateral stability derivatives of a linear aircraft. The use 
of four different types of control deflections as forcing functions pro- 
vided a means for matching of computer responses which corresponded 


to a wide range of flight maneuvers. 
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The method is considered suitable for further investigation, and 
development for use under different flight conditions and with non- 


linear aircraft. 
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INTRODUCTION 


Since the close of World War I], and the advent of the high speed 
and altitude characteristics of jet aircraft, the need for a quantitative 
knowledge of the dynamics of flight has resulted in rapid development 
of the field of aircraft dynamics. The past six or seven years, in 
particular, have brought considerable progress in the field of dynamic 
stability flight testing of aircraft and missiles. The early theoretical 
and experimental investigations in this area by the Cornell Aeronau- 
tical Laboratory, and by the Instrumentation Laboratory of the Massa- 
chusetts Institute of Technology, have been supplemented on many 
fronts, and in particular, by the NACA, which in recent years has 
developed its dynamic stability testing techniques to the point where 
free-flight testing of rocket powered and telemeter equipped models 
has become commonplace. 

The role of full scale dynamic flight testing is twofold. First, 
it provides a designer with the conventional data on the characteristics 
of an aircraft for the particular conditions of the test. Secondly, and 
in the long view, more importantly, dynamic flight testing provides a 
means for solution of the inverse problem of dynamics, wherein 
the equations describing a system are assumed known, its output in 
response to suitable forcing functions is measured, and thus the co- 


efficients of the equations, or stability derivatives, can be determined. 





Once established, these stability derivatives are of considerable value. 
They are indispensable to the aerodynamicist who must calculate sta- 
bility characteristics over a wide range of flight conditions, and who is 
thus provided with a measure of the accuracy of the theoretical and 
empirical derivatives used in preliminary design calculations. Deter- 
mination of the stability derivatives from flight tests also aids the 
designer by indicating lines along which the stability may be improved, 
either by geometric design changes, or by artificial augmentation. In 
the latter regard, the successful design of any system for automatic 
control of an aircraft is greatly dependent upon an accurate knowledge 
of the stability derivatives of the airplane. 

The basic approaches to the problem of the determination of 
stability derivatives from dynamic flight testing, and in particular, 
from experimental time histories, can be separated into two categories 
- the "equations of motion" method, which essentially solves the in- 
verse problem described above, anda second category, "response 
curve fitting," which identifies a collection of methods using various 
mathematical and graphical concepts to extract from the data certain 
combinations of stability ee batives and inertia terms called transfer 
function coefficients. All of the individual methods which make up the 
above two categories will not be described in this report, but an excel- 
lent discussion of each current method can be found in Ref. 1. Table I 


presents a summary and comparison of these various methods for 





analysis of dynamic systems. As can be seen from Table I, the choice 
of a method for analyzing dynamic response data will depend primarily 
on whether transfer functions, or actual stability derivatives are 
desired. 

Several investigations have been reported (ie, Ref. 2) which 
have attempted to obtain stability derivatives by applying a particular 
combination of the above referenced methods which essentially involved 
the conversion of time response to frequency response, followed by a 
matching process which utilized both dynamic response templates and 
an analog computer to produce transfer function coefficients. These 
constants were then shown to uniquely determine the time response of 
an analog of the test airplane which accurately duplicated the time his- 
tory of the actual test vehicle. Evaluation of the actual aerodynamic 
stability derivatives, however, required supplemental calculations, and 
was possible only in the case of those derivatives which exerted the 
strongest influence on the airplane response. 

The major limitations of the above method are the length of time 
required for conversion of the flight data from time response to fre- 
quency response form, and the necessity for the subsequent use of tem- 
plates to determine approximate transfer function coefficients. It was 
considered that both of these limitations could be avoided in the applica - 
tion of the above method if the matching process were to be applied 


directly to the transient flight data. 
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This report, then is concerned with a study of the feasibility of 
obtaining certain lateral stability derivatives by matching the transient 
responses from an analog computer on which the operational equations 
for a NAVION aircraft were set up, to the actual airplane lateral re- 
sponse in sideslip angle, roll rate, and yawrate. The matching pro- 
cess was achieved by supplying time histories of the actual recorded 
control motions as inputs to the computer, and systematically varying 
those coefficients of the analog which corresponded to the stability 
derivatives of the dynamic equations of motion. 

An additional objective of this investigation was to determine the 
suitability of a pulse width telemeter system as a data collecting 
scheme for a dynamic flight test program requiring very accurate 
determination of the recorded quantities. 

This investigation was carried out at the James Forrestal Research 
Center, Princeton University, Princeton, New Jersey, during the 1956-57 


academic year. 





EQUIPMENT AND PROCEDURE 


EQUIPMENT 
Airborne Instrumentation 

Test Airplane. The test airplane was a North American NAVION, 
which is an all metal, four place, low wing, single engine airplane. 

The NAVION is equipped with a Continental E-185 engine, rated at 185 
HP at 2300 RPM and 29 inches manifold pressure, at sea level. The 
wing is unswept, and the tail configuration conventional. The control 
surfaces are all of beaded skin construction, and have no trailing edge 
extensions. Aileron and rudder trim tabs are of the fixed bend type. 
The elevator tab is of conventional design, and can be adjusted from the 
cockpit. 

The test airplane was modified by the removal of the two rear 
cockpit seats to provide space for the test instrumentation. (Fig. 6). 
The test airplane had a take-off gross weight of 2870 pounds with full 
service, pilot and co-pilot, and all test equipment. 

A three view drawing and photograph of the NAVION are shown 
as Fig. 1 and Fig. Z, andthe general specifications of the test airplane 


are tabulated in Table 7. 


Power Supply. Power available for instrumentation in the test 
airplane included both 12 and 28 volt direct current, supplied by a 28 


volt engine driven generator, andtwo 12 volt batteries connected in 











series. A 500 volt-ampere inverter, with both frequency and voltage 
regulation, was installed behind the co-pilot's seat to provide 115 volt 
400 cycle three phase alternating current power to the rate gyros and 
the telemeter transmitter. 

Considerable effort was expended in obtaining an inverter which 
could be frequency controlled, since the accuracy of the rate gyros was 
dependent upon the input frequency being the same as the calibration 
frequency of 400 cycles per second. 

A circuit diagram of the aircraft 28 volt system and the inverter 


system is included in Fig. 3. 


Aircraft Instrumentation. Standard aircraft instruments were 
used to determine indicated airspeed and altitude. Measurement of the 
following quantities was required with respect to stability axes: side- 
slip angle, yaw rate, roll rate, and aileron and rudder position angles. 
In choosing the type of components needed to measure the above quan- 
tities, an attempt was made to utilize instrumentation compatible with 
the telemetering equipment. For this reason, the transducers selected 
were of the potentiometer type, not only because they were capable of 
providing signals of a sufficient magnitude, but also because their excita- 
tion source readily provided zero and full scale reference voltages. 

In view of some past history in this type of flight testing of various 
amounts of coupling between the response of the aircraft and certain of 


the instrumentation components, considerable attention was given to the 





choice of transducers with satisfactory frequency and phase response 
characteristics. These remarks particularly apply to the choice ofa 
suitable device for recording sideslip angle. 

In this investigation, sideslip angle was sensed by a light vane 
mounted on a boom approximately one maximum fuselage diameter 
(four feet) ahead of the left wing tip. The moment of inertia of the vane 
was very low, and no gearing was introduced in the linkage in order to 
keep the natural frequency and damping within acceptable limits. An 
Electro-Mechanical Laboratory ultra-low torque potentiometer, of only 
.008 inch-ounces friction torque, was obtained for use with this side- 
slip vane. In Appendix A is presented a detailed analysis of both the 
theoretical and experimentally obtained natural frequency and damping 
ratio of the sideslip vane. 

The sideslip vane had originally been designed to operate through 
an autosyn transmitter, and thus required some modification in order to 
accommodate a potentiometer. 

In order to reduce boom resonance, as had been experience in 
previous tests with the NAVION, the boom for the sideslip vane was 
stiffened by three angle strips equally spaced around the periphery, and 
running the full length of the boom. The sideslip vane and boom are 
shown in the assembled condition in Fig. 4. 

Rudder and aileron angles were measured by 10,000 ohm Model 


5301 single turn Helipot precision potentiometers. The rudder pot 





was attached to the rudder itself, and was actuated by a cable under 
spring tension. At the time the instrumentation system was designed, 
the flight test program was planned to include control inputs from the 
rudder only. For this reason, the instrumentation installed for re- 
cording aileron angle was intended only to indicate if the ailerons 
remained locked throughout the rudder input and the resulting tran- 
sient oscillation. Accordingly, the pickoff of aileron angle was made 
from control] cable movement, with the pot attached close to the right 
aileron in order to avoid the effects of cable stretch. As is discussed 
in Appendix C, the later decision to include aileron deflections as 
forcing functions required subsequent correction of the recorded aileron 
data. 

The measurement of yaw and roll rates was made by type 36128 
Giannini rate gyros. This type of gyro operates on a power supply of 
115 volts 400 cycle three phase AC to the motor, and 115 volt single 
phase AC tothe heater. The transducer is of the potentiometer type, 
with a DC resistance of 5250 ohms. The gyros were installed as near 
the aircraft center of gravity as possible. However, since the gyros 
were of the angular velocity rather than acceleration type, there was 
not the necessity for precise location that is associated with the latter 
type gyro. A level and inclinometer were used to align the gyros so 
that the input axes would be parallel to the appropriate body axis of the 


aircraft. Since deflection of the gyro spring-restrained gimbal was 





only two degrees at full range, cross-coupling effects were considered 
to be negligible. However, in order to minimize any error due to rate 
coupling, the spin axis of each gyro was aligned parallel to the body 
Y-axis. The gyros were bolted directly to a rigid fuselage member. 
Since the engine vibration frequency was sufficiently above the natural 
frequency of the gyros (17 cycles per second), there was no engine 
noise introduced into the gyro signal. These gyros were obtained on 
loan from the Instrumentation Laboratory of the U. S. Naval Air Test 


Center, NAS Patuxent River, Maryland. 


Signal Circuit Design. The component of the instrumentation 


which controlled the design of the signal circuit was the sideslip vane. 
The requirements for recording the sideslip angle without introducing 
phase differences were negligible torque, anda relatively high natural 
frequency of the vane. The transducer selected has been shown to have 
met these specifications (Appendix A), but the fact that it was a 360 
degree single turn potentiometer pickup, whereas the sideslip vane was 
expected to operate through only a moderate arc, necessitated the use 
of 45 volt DC excitation of the sideslip pot rather than the zero to five 
volts range supplied to the remaining transducers. Since the resolution 
of the particular sideslip pot used was exceptionally good, the accuracy 
of the sideslip system was comparable to that of a system utilizing 
autosyn pickups, but had none of the problems of AC signal circuits 


and was completely compatible with the telemetering and recording 
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equipment. The use of a 45 volt DC signal circuit supply, however, 
necessitated a method of supplying this voltage to the sideslip circuit, 
and at the same time, limiting the full scale voltage in the remaining 
Signal circuits to five volts. 

In achieving the above, it was desirable to utilize a common DC 
voltage source for all transducer pickups, since the telemeter was to 
sample the particular supply voltage available 20 times per second, 
and correct this to full scale by a continuous calibration. 

In addition to the above requirements on the signal circuit design, 
it was found necessary to bias the sideslip voltage pot in order that the 
sideslip signal voltage would also vary between zero and a full scale 
value of approximately five volts. 

The particular signal circuit design employed is shown in Fig. 5, 
in schematic form. The essential element of the circuit is the bias 
supply voltage potentiometer, pot "A". With the sideslip vane against 
the limit stop in one direction, pot "A" was adjusted so that the output 
voltage of the sideslip pot was zero. The vane was then deflected to the 
opposite limit stop, and the output sideslip pot voltage recorded. Pot 
"B*" was then adjusted to bring the full scale reference voltage, and 
thus also the voltage in the remaining signal circuits, in agreement with 
the full scale sideslip voltage. This procedure was repeated two or 
three times, as necessary, in order to nullify any loading on the system. 


Once adjusted, however, pots "A" and "B" were locked, and required 


no further adjustment for the remainder of the project. 





During the initial flight and transmission tests, considerable AC 
ripple was found to enter the DC signal voltage circuits, very likely 
from the 400 cycle inverter. Condensers C and D in Fig. 5 
were installed, and were found to reduce the ripple to a negligible value 
without otherwise altering the form or amplitude of the signals. 

Fig. 6 shows the general location and arrangement of the air- 


borne instrumentation components. 


Telemeter System. The telemeter unit was of the pulse width 
type, and was designed and built by the Applied Science Corporation of 
Princeton (ASCOP). This unit was the key element in the instrumenta- 
of this project. The basic functional components of this pulse width 
system are transducers, all of the potentiometer type; a time division 
multiplexer, or commutator; and electronic coder, or keyer to put the 
intelligence into pulse form; an RF transmitter; and ground station 
equipment for receiving, decoding, and recording the data in usahle 
form. 

The airborne component of the system is designated the ASCOP 
D Series PW Multicoder, and the ground station equipment the ASCOP 
M Series PW Ground Station. A schematic diagram of the functions 
of the pulse width data system, as adapted for use with this particular 
project, is presented in Fig. 7. 

A unique feature and very great advantage of pulse width data sys- 


tems is that due to the multiplexing or time sharing operation, all data 
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channels traverse the same path during transmission. Since zero and 
full scale references are transmitted with the data, this results in a 
calibration signal, effective for all channels, being transmitted 20 
times each second. These signals supply a continuous two-point cali- 
bration, and thus allow for the removal of almost all of the system 
errors, with the exception of errors introduced in the instrument pick- 
ups and the mechanical switch, and they make possible the use of this 
telemeter in tests of long duration and varying environmental conditions. 
The output voltages from the aircraft instrumentation pickups, 
and the calibration voltages described above, are sampled in sequence 
by a rotary switch. At the normal sampling rate of 20 samples per 
second on each of the 45 channels, the data pulse length for any one 
channel is short in comparison to the interval from sample to sample. 
This permits the use of a storage circuit to create a nearly continuous 
voltage output. The output of the sampling switch is a series of data 
pulses whose voltage amplitude is proportional to the amplitude of the 
physical quantity measured by the corresponding pickup. These samples 
are supplied to an electronic coder circuit (keyer) which converts the 
amplitude pulses into pulses having time durations proportional to the 
measured quantities. These pulses are used to produce modulation of 
the transmitted RF carrier. The RF signal is detected at the ground 
station, with the output of the receiver being a train of pulses identical 


in nature to the output of the converter in the airborne equipment. 
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These pulses are separated by a switch running synchronously with 

the switch at the transmitter, with the pulses representing individual 
data channels being placed on separate outputs. The pulses are then 
demodulated to yield continuous voltages representing the intelligence 


transmitted. 


Ground Station Instrumentation 

The ground station equipment included the Telemeter Ground 
Station (described above), a Magnetic Tape Recorder, a GEDA 
Analog Computer, and suitable recorders. These components are 
illustrated in Fig. 8. 

Magnetic Tape Recorder. An Ampex Model 309C dual track 
three speed tape recorder was an integral component of the ground 
station equipment, and provided a convenient means for recording, 
storage, and playback of the standard pulse width modulated data. 
Provisions were also included for placing voice transmissions directly 
on the magnetic tape from either or both the aircraft and the ground 
station. 

GEDA Analog Computer. The analog computer used in this 
investigation was a Goodyear Aircraft Corporation Model L3 (GEDA) 
linear electronic differential analyzer. Twenty-four automatically 
stabilized DC computing amplifiers were available with open-loop DC 


gain of greater than 5 x 10°, and negligible drift. The computer 
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incorporated an automatic error indicator, and had a guaranteed 
accuracy of one percent. Provisions were available for accurately 
setting computer board potentiometers by the use of a special null 
indicator. 

Data Recorder. Qualitative analysis of the flight data during 
tape playbacks was made ona Brush two channel recorder. Once the 
particular data run was selected for further analysis, a Sanborn Model 
154-100B four channel recorder was used. This instrument had a very 
high natural frequency (42 cps), and the low frequency response was 


flat to zero cps. 


Instrument Calibration 

Five components of the instrumentation system required calibra- 
tion. These were the two rate gyros, the two control deflection poten- 
tiometers, and the sideslip vane potentiometer. The aircraft altitude 
and airspeed indicating systems were not calibrated, as previous cali- 
brations indicated that the errors in these systems were negligible for 
the conditions of this investigation. 

The two rate gyros were calibrated on a pendulum which consisted 
of a 3/4 inch diameter aluminum shaft and a 12 inch diameter instru- 
ment pan. The pendulum, the overall length of which was 121.75 inches, 
was pinned to the overhead truss structure of the Forrestal hangar. 


A scale was mounted on a small table in such a way that displacement 
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of the pan from its rest position could be measured. The gyros were 
mounted one at a time on a wooden base plate which was clamped to 

the pendulum pan, and their potentiometer outputs were connected to 

a Brush Pen Recorder. The pendulum was displaced from its rest 
position and allowed to swing freely. Because of friction in the system, 
swing amplitudes were measured on two successive swings in the same 
direction and averaged. The two maximum voltages recorded during 
the swing were also averaged and the maximum rate attained during 


the swing was determined from the following relation: 


: Re ao 
OF. - BI = 
where: 


P = natural period of the pendulum, determined from the 
Brush Recorder trace, seconds. 


Bi = length of the pendulum, inches. 

d <= average maximum displacement of the pan from rest 

position, inches. 

The pendulum was allowed to swing until once again at rest with 
readings taken as described above over the entire range of gyro rates. 
The angular rate gradient of both gyros was found to be linear and equal 
to 0.610 degrees per second, per percent full scale potentiometer 
voltage. The gyro calibration curves are presented as Fig. 9. It was 


found that the gyros were very sensitive to the frequency of the 400 





16 


cycle supply voltage and as a result, great care was taken during 
calibration, and also during the actual flight test program, to main- 
tain the frequency of the supply voltage at exactly 400 cycles per 
second, 

The control deflection potentiometers were calibrated by meas- 
uring the resistance between the output terminals over the full range 
of control deflections. The actual control deflections were measured 
with a vernier protractor capable of accuracy to 0.1 degree. The 
rudder deflection gradient was found to be linear and equal to 0.405 
degrees per percent full scale potentiometer voltage. The aileron pot 
calibration was based on total aileron angle, and was performed over 
the entire range of aileron travel in order He Gani the extent of 
aileron differential for the NAVION aircraft. The calibration was 
then converted to average aileron angle, and was found to be linear at 
a constant value of 0.1933 degrees per percent full scale, throughout 
the deflection range used in the flight test, which included -10 to 10 
degrees. The amount of aileron differential gearing in this range was 
found to be extremely small. 

The sideslip vane potentiometer was calibrated in the same 
manner as the control deflection potentiometers, and the calibration 
curve is presentedin Fig. 10. The sideslip gradient was found to be 
linear and equal to 0.468 degrees per persent full scale potentiometer 


voltage. 
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Since the telemeter equipment corrected its output to percent 
full scale voltage and since full scale voltage on all the signal circuit 
potentiometers varied by the same increment as the reference voltage 
for small fluctuations of the reference, it was unnecessary to cali- 
brate the potentiometers through the telemeter equipment. 

Careful consideration was given to the possible errors which 
might be introduced in the data by the dynamics of the recording in- 
struments. The natural frequency of the rate gyros was 17 cycles 
per second and their damping ratio was between .45 and .7. The 
natural frequency and damping ratio for the sideslip vane were experi- 
mentally determined in flight to be 15.44 cps and .081 respectively. 
(See Appendix A.) Since the natural frequency of the aircraft was 
approximately 0.4 cycles per second, the attenuation of the signals 


due to the dynamics of the instruments was effectively zero. 
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PROCEDURE 


General 

Aircraft transient responses in the lateral modes of motion were 
obtained by introducing rudder and aileron control motions of the pulse, 
doublet, and step-function types. Rudder pulses, rudder doublets, and 
aileron doublets were found to be the most effective in exciting the 
® Dutch Roll" motion, which in turn produced the most satisfactory 
transient time histories for the investigation in question. Aileron step- 
functions were also used, however, primarily to provide the time 
history of roll-rate response to lateral control motion. 

Ideally, sufficient transient response data for all types of forcing 
functions could be obtained in one flight of approximately an hour's 
duration. However, due primarily to the need of the investigators to 
gain some proficiency with the instrumentation and technique of applying 
satisfactory forcing functions, several flights were made in the process 
of collecting data. It is worthy of note that the basic scheme of instru- 
mentation and data collection proved singularly reliable during the flight 
test - data collection phase. After initial installation the instrumentation 
required only a very minimum of modification and adjustment, and in 


each instance the trouble was minor in nature and easily remedied. 
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Aircraft 

Preflight preparation consisted of checking the full scale voltage 
applied to the sideslip potentiometer to insure that it matched that of 
the other signal potentiometers, and to insure that the full scale signal 
voltage did not exceed the 4 - 6 volt range which the telemetering 
equipment needed for automatic scale calibration. The fuel tank was 
checked and if necessary topped-off to insure that the aircraft weight 
at the beginning of each flight was constant. When deemed advisable, 

a complete ground check of all communication and signal transmission 
circuits was made prior to take-off. After completing the ground 
checks all test equipment and the aircraft generator were secured to 
insure maximum engine power for take-off. 

During the climb to test altitude, the generator and inverter were 
energized and power was applied to the rate gyro heaters and to the air- 
borne telemeter heater and filament circuits. Upon reaching test alti- 
tude, voice contact was established with the telemeter ground station 
and all instrumentation circuits were energized. Because of an appar- 
ent deficiency in the inverter frequency control, the inverter frequency 
was monitored continuously and maintained at 400 cps by varying the 
DC voltage supplied by the aircraft generator to the inverter. Voltage 
control was obtained by means of an instrument-panel mounted potenti- 
ometer which was tied into the aircraft voltage regulator circuit (Oper- 


ation of the synchronous motors in the rate gyros at exactly 400 cps was 
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necessary to insure that the output signals followed the calibration 
curves obtained at that frequency). 

All flight tests were conducted at a density altitude of 5000 feet 
and at a true airspeed of 129.2 mph (120 mph IAS at 5000 density 
altitude). A large plot of free air temperature versus pressure alti- 
tude was carried in the cockpit to facilitate adjusting pressure altitude 
(altimeter set at 29.92 in. Hg) with temperature to obtain and main- 
tain 5000 feet density altitude. 

The general technique employed in the actual test runs fits the 
following description. A condition of steady-state flight was established 
at test altitude and airspeed and held for at least 30 seconds prior to 
each control input. This period of steady-state flight insured that the 
aircraft was as nearly trimmed and established in steady state condi- 
tions as possible, and provided a sufficient period of recorded steady 
state signals to allow for circuit zero-signal biasing when the data was 
played back into the analog computer or pen recorder. At the end of 
the steady state period a radio transmission identifying the type of con- 
trol displacement to be made was sent and recorded on the magnetic 
tape. Simultaneously, all controls were locked mechanically in the 
cockpit with the exception of the control to be displaced. Five seconds 
after the above transmission the control forcing function was applied 
and the control was then returned to the original steady-state position 


and locked. All controls were held locked until the resulting aircraft 
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oscillation had essentially subsided. In the case of the aileron step 
functions the displacement was applied as rapidly as possible and the 
control was then held as rigidly as possible. The resulting motion 
in this case was permitted to proceed for 4 to 6 seconds before 
recovery was initiated. 

To lock the elevator and ailerons when forcing the aircraft with 
rudder displacements, a chain and turnbuckle system attached to the 
deck on the right side of the cockpit was used. The copilot's control 
wheel was removed and replaced by a sprocket wheel with a set screw 
to lock it to the control shaft. When the aircraft was established in 
steady-state flight, the sprocket wheel was pushed forward on the con- 
trol shaft until it butted against the instrument panel. The chain was 
then put on the sprocket and the turnbuckles taken up, making the chain 
taut. Just prior to initiating the rudder deflection, the set screw was 
tightened, locking the sprocket fitting to the control shaft. With the 
chain on and tight, the ailerons were effectively locked, and, with the 
sprocket fitting against the instrument panel, no forward motion of 
the control shaft was possible. The copilot then held forward pressure 
on the sprocket fitting to prevent aft movement of the control shaft. 

To lock the rudder when putting in aileron displacements, and to 
return the rudder to its steady-state position after rudder displace- 
ments, a chain and turnbuckle system on the pilot's right rudder pedal 


was employed. This was attached to the right hand seat track on the 
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pilot's side. Because the aircraft required a slight amount of left 
rudder for trim at the test flight condition, this attachment was also 


useful in establishing trimmed flight by adjusting the turnbuckle. 


Ground Station 

In order to visually monitor the data runs, an external patch- 
board provided on the telemeter receiver (See Fig. 8) was wired in 
such a way as to utilize a voltmeter-type translator for each of the 
recorded variables, and also to allow the decoded output to be placed 
on apen recorder. The telemeter receiver and the magnetic tape 
recorder had been integrated prior to this investigation, thus requiring 
only the proper positioning of controls on the telemeter console to 
accomplish tape recording of the flight data. Each quantity measured 
during the flight tests was transmitted on three separate telemeter 
channels. Thus, each measured variable was sampled 60 times each 
second. Once the airborne equipment was operating, the telemeter 
receiver was tuned for maximum signal reception, and the "zero and 
"full-scale" reference voltages were adjusted for proper calibration of 
output signals. 

All pertinent communications between the ground station and the 
aircraft were recorded on magnetic tape along with the flight test data. 
Prior to commencing the recording of test data on each flight, the 
ground station operator pre-recorded on the magnetic tape the data 


necessary to identify the particular flight, and test conditions as 
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regarded aircraft weight, lift coefficient, and the telemeter channels 
being used for each measured quantity. 

Concurrent with the reception and tape recording of the flight 
test data, a two or four channel pen recorder was used to present 
visually for monitoring all or a part of the test data as the maneuvers 
were being performed. Thus, if the control input did not produce well 
excited responses, or if a control] surface that was to be held locked 
moved inadvertently, or if the forcing control was not returned to, and 
held at, the original steady-state position this was immediately appar- 
ent. In any of these eventualities the run could be rejected immediately 


on a qualitative basis and the test repeated. 


Data Reduction 

The only flight variable requiring correction to the measured 
quantity was sideslip angle ( & ) which had to be corrected for location 
of the sideslip vane ahead of the aircraft center of gravity. The required 


correction was: 


Pome. = J xmas a x (y) 


In order to have all flight variable signals in the same angle- 
time-voltage relationship, the flight data was processed through scaling 
circuits on the analog computer prior to plotting on a pen recorder. 

The scaling factors used were determined from the calibration curves 


for each of the signal potentiometers. Where necessary, the signs of 
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the measured variables were changed to correspond to the sign of the 
corresponding output of the analog probelm. A .25 uf capacitor was 
placed across each of the output terminals of the telemeter. This 
size capacitor was found to help in eliminating high frequency noise in 
the flight data signals while at the same time having negligible effect 
on the amplitude and phase of the signals. The computer circuits 
used in correcting the sideslip angle and scaling the flight test data 


were as follows: 
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Flight Data Scaling and Correction Circuits 


In order to use the actual flight test control deflections as forcing 
functions for the analog problem it was necessary to feed them into the 
computer through a bias circuit so that the steady-state control position 


signals, which were of the order of -50 volts at the telemeter output 
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terminals, could be reduced toa zero voltage input to the analog 
problem. A potentiometer was also required in each input circuit to 
adjust the amplitude of the control deflection voltage to match the 
scaled magnitude of the control deflection on the flight test data plots. 
As in the data reduction circuits, a .25 uf capacitor was put across the 
telemeter output terminals. The above described circuits are included 


in Fig. 13. 
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DISCUSSION 


General Description of Overall Program 


In obtaining flight data, four different types of control deflections 
were used to excite the lateral modes of the NAVION. These were a 
rudder pulse, a rudder doublet, an aileron doublet, and an aileron step 
function. By means of telemeter equipment, the transient response of 
the aircraft as well as the control deflections as a function of time were 
recorded on magnetic tape at the ground station. From a large selec- 
tion of such data, one representative trace was made of each of the four 
types of forcing functions mentioned above. This was accomplished by 
processing the flight data through scaling circuits in the analog computer 
and recording the variables in time history form on a Sanborn recorder. 
A transparent overlay was then made of each of the four traces of 
scaled flight data in order to facilitate comparison of the flight and com- 
puter responses. 

Next, the lateral stability derivatives of the NAVION were 
theoreticallyestimated. The lateral equations of motion based on these 
derivatives were placed on an analog computer in such a way as to enable 
independent variation of each stability derivative. The magnetically 
recorded signal corresponding to one of the four chosen control deflec- 
tions was then fed directly into the computer to force the analog. Both 


the character and magnitude of the forcing function signal were checked 
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by means of the overlay of the flight data. Once the forcing function 
was properly adjusted, the computer responses to the forcing function 
were recorded and compared to the flight responses by means of the 
same overlay. Adjustments were then made to the computer deriva- 
tives as necessary to improve the match between the flight data and 

the computer response. The procedure of forcing the computer was 
then repeated, using another one of the four forcing signals, and mak- 
ing appropriate adjustments to derivative values. This procedure 
proved to be a converging iteration process when repeated several 
times, using all four types of forcing functions. Whena set of deriva- 
tives was found such that the response of the computer gave satisfactory 
matching with the flight data for all four of the forcing functions, and no 
further improvement seemed justifiable in relation to the amount of time 
involved, the value of each stability derivative was read out of the com- 


puter and tabulated. 


Preliminary Considerations 
Theoretical Background. At the onset, of the investigation, it 


was necessary to make a purely theoretical analysis of the proposed 
method of attack of the problem. In effect, the pertinent questions were: 
a) If a set of coefficients is determined which, when placed in the 
equations of motion, duplicate the response of the aircraft, 
can it be proved that these coefficients are unique rather than 


just one set of an entire family of possible solutions? 
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b) With two possible methods of exciting the lateral mode, i.e., 
the rudder and the ailerons, how much data should be col- 
lected and what should be the character of the control deflec- 


tions used to excite the motion? 


Examining these questions, it was determined that one was related 
to the other. The answer to question (a) is a function of the amount of 
independent data available. Should there be a case where, in fact, an 
entire family of solutions exists, then there would be an interdependence 
between at least two of the coefficients, and one coefficient could be ex- 
pressed in terms of one or more of the other coefficients. If it were 
then possible to obtain more independent data, a point would be reached 
where there were sufficient conditions to limit the family of solutions to 
a Single solution. Applying the latter case to the investigation at hand, 
it would then be known that if a set of coefficients were obtained by the 
response-matching technique, these coefficients would be unique, i.e., 
the only set of coefficients which would produce the matching of re- 
sponses. The stability derivatives contained in the lumped coefficients 
would then indeed be the stability derivatives of the aircraft. The 
amount of independent data available would be in part a function of the 
control deflections used, and in addition would depend on the number of 
quantities such as roll rate, yaw rate, and sideslip angle which were 
recorded during the response. 


Fortunately, a theoretical investigation into this subject was 
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available in Ref. 4. The authors, Messrs. Klein and Sedney, assume 
that the equations of motion are linear with constant coefficients (how- 
ever, it is pointed out that the same method of analysis can be applied 
to more general types of equations), and also that the data is free from 
experimental error. The results of their analysis is a simple working 
rule which can be stated in words as follows: In each equation describ- 
ing a system, the number of unknown coefficients which can be deter - 
mined plus the number of known coefficients which may be checked are 
equal to the total number of measured quantities less the number of 
equations. 

This result was applied to the lateral equations of motion of the 


NAVION, with J,, Jz, J 


xg» and C, known, The quantities which 


were to be recorded during the testing were roll rate, yaw rate, side- 
slip angle, and control deflections. In application of the Klein and 
Sedney criterion, roll angle, roll acceleration, yaw angle, yaw accel- 
eration, and sideslip rate are considered to be recorded in a concealed 
form in the other time histories, since they are obtainable by differen- 
tiation or integration of the appropriate trace. By following the working 
rule of the criterion, it was found that in each of the three equations of 
motion there was actually an excess of independent data under the re- 


cording conditions proposed, and therefore the solution for the coeffi- 


cients would be unique. 
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Practical Aspects. Although no theoretical limitation exists in 
the case of the lateral mode, there is a practical limitation which 
applies to the analysis of all dynamic response data. The following 
statements are extracted from Chapter 10 of Ref. 1: 


Instrument errors and unknown atmospheric disturbances 
are usually present to a sufficient degree to introduce significant 
errors into the solution of the inverse problem. As a general 
qualitative rule, the accuracy with which one can obtain each of 
the stability derivatives is associated with the importance of that 
stability derivative in influencing the motion of the aircraft dur- 
ing the test, the highest accuracy being obtained for the more 
important stability derivatives. 


Consideration of the means for improving the accuracy 
with which the less important derivatives can be obtained leads 
to these three approaches to the problem: 

1. Overall improvement of the accuracy of the recording 
instrumentation and improvement of the calibration 
and pre-flight testing techniques. 

2. Selection of test maneuvers which increase the force 
and moment contributions due to the smaller deriva- 
tives compared to the forces and moments due to the 


other derivatives. 


3. Performance of flight tests under so-called ideal 
conditions, i.e., air is smooth and stable. 


Selection of flight test instrumentation components was made 
with the requirements of (1) above inmind. As previously discussed 
under “Airborne Instrumentation, the transducers were chosen to be 
compatible with the recording instrumentation in order to eliminate 
excessive electronic modification of the signal voltages prior to record- 


ing. This effectively reduced possible sources of error. Calibration 
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was carefully performed and the resolution and linearity of the com- 
ponents are readily apparent from the calibration curves, as presented 
in Fig. 9 and Fig. 10. Test flights were performed under conditions 
of smooth and stable air. 

The selection of suitable forcing functions for the lateral mode 
was initially considered from the viewpoint of Klein and Sedney, who 
state that theoretically, the response of the aircraft to a single forcing 
function is sufficient to establish the stability derivatives, with the 
exception of the control derivatives. Accordingly, a rudder pulse was 
selected as the forcing function most suitable for excitation of the 
lateral mode. However, the prospect of extending the method of analy- 
sis to obtain all of the control derivatives, plus the possibility of 
improving the accuracy of certain derivatives through the use of addi- 
tional forcing functions, was made very attractive when considered 
from the standpoint of paragraph (2) of the above quotation from Ref. 1. 
Four forcing functions were therefore decided upon. Considering the 
rudder pulse as the basic data, a rudder doublet was utilized to in- 
crease the amplitudes of the lateral responses, and to reduce the 
tendency of the aircraft to fall off on one wing during the transient as 
had been observed in the case of excitation by a rudder pulse. The 
aileron doublet was considered as an effective means of evaluating 
dh , Cn, , and the aileron derivatives, whereas a time-history of 
the roll rate response to an aileron step function was considered to be 


the best means for analyzing the stability of the spiral mode. 
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Analog of the Test Aircraft 
Theoretical Approximation of Derivatives. In order to construct 


an analog of the lateral equations of motion, it was necessary to esti- 
mate the value of the stability derivatives from whence the approximate 
value of each lumped coefficient could be determined. Anticipating the 
possibility of extremely slow convergence of the iteration process used 
on the computer, the theoretical value for each derivative was deter- 
mined as accurately as possible by several different methods and then 
a median value was taken. The results and a complete explanation of 


this work is contained in Appendix B. 


Lateral Equations of Motion. The lateral equations of motion of 


a linear aircraft are as follows: 
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For ease of construction of the analogy the equations can be 
rewritten incorporating the following changes: 

a) in terms of the highest derivative. 

b) in real time by the explicit insertion of the non-dimensional 


time constant, 2. 
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c) for simplicity, the algebraic signs changed where necessary 


so that all work from this point on could be conducted in 


terms of the absolute value of each derivative. 


The resulting equations are: 
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The numerical values of the non-dimensional constants and the 
estimated derivatives were then substituted into the equations. The 


terms containing J,, were found to be so small in comparison to the 


other terms that they were immediately discarded as negligible: 
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The constructed analog for the above equations is shown 


schematically in Fig. ll. 
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Development of Method of Analysis 
Prior to the analysis of the actual flight test data, a step by 


step procedure was developed for matching the computer response to 
flight data. First, an investigation was made to determine the effect 
of each stability derivative on the response of the system in terms of 
amplitude, period, and damping. The analog of the test airplane 
(Fig. 11) was placed on the analog computer, and the potentiometers 
set to correspond to the NAVION theoretical stability derivatives of 
Table V. An electronic function generator was used to force the 
analog with a square pulse type function, to approximate the form of a 
rudder pulse. A series of runs was then made during which the sta- 
bility derivatives were varied successively by a factor of two and then 
returned to their original theoretical values. From these tests, the 
change in the response due to the change of any one derivative was re- 
corded, and the results included in Table III-a. A cross-table was 
then prepared showing which derivatives should be varied to obtain a 
desired change in the response of the system. These results are 
shown in Table III-b. 

With the aid of the information presented in Tables III-a and 
Iil-b, an iterative procedure for the matching of computer response with 
flight data was outlined as follows: 

(a) Compare flight data and computer response curves obtained 


from the same forcing function. Note discrepancies in 





(b) 


(c) 


(d) 


(e) 


(f) 


amplitude, phase, damping, etc., between corresponding 
responses. 

Concentrate on the major discrepancy. By reference to 
Tables IJIJ-a and I0-b, determine which derivative, or 
derivatives, could be altered to correct this type of 
discrepancy. 

Iteratively apply the corrections indicated in (b) above, 
commencing with the derivative having the most pronounced 
effect on the particular discrepancy. When no further 
improvement can be observed by adjusting the value of the 
first derivative, continue in a consecutive adjustment of all | 
the derivatives affecting the discrepancy. 

Repeat (a), (b), and (c) above for additional discrepancies 
until general agreement is obtained between flight data and 
the computer response curves. 

Continue by comparing flight data and computer response 
curves obtained from a second forcing function. Again, 
repeat (a), (b), and (c) until agreement is obtained. 

Repeat (e) for as many forcing functions as are available or 
practical. When the response to all forcing functions have 
been examined once, return to the first forcing function and 
repeat the entire process outlined above; continue until a 
set of derivatives is obtained which give the best possible 


response match to all forcing functions used. 
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Preliminary Test of Feasibility of Method 


A program was next performed to establish the feasibility of 
using an analog computer to curve-fit flight data and yield unique 
stability derivatives. It was considered that if the transient response 
of the computer to a rudder-pulse type of forcing function were record- 
ed, and considered as error-free ‘flight data, and the computer 
potentiometers corresponding to stability derivatives then changed, 
that subsequent re-matching of the original synthetic "flight data - 
response should yield computer potentiometer settings identical to 
those of the - flight data. It was found that this type of analysis 
produced excellent results, with the exception of the determination of 
Ch. ; Cy, » and Cys - The responses were found to be too 
insensitive to changes in these three derivatives in order to establish 
them with satisfactory accuracy. The final value of C nb contained 
an error of approximately ten percent. All other derivatives were 


obtained within an error of five percent. C and Oh, were 


8 
obtained with virtually no error; this latter result was anticipated, 
from consideration of the large relative magnitude of the lumped 
coefficients containing these terms. 

Unfortunately, the above program was not extended to the use of 
an aileron input. Although the results were gratifying for the rudder 


forcing function used, they might well have been improved if the addi- 


tional forcing function data had been utilized. In particular, the use of 
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aileron forcing functions might have resulted in a better indication of 
the effect of the derivative ys : 

The investigators considered that the results of this preliminary 
program established the feasibility of the basic method. It was also 
encouraging to note that with a minimum of practice, the iteration pro- 
cedure developed for matching the analog responses to flight data 
resulted in moderately rapid convergence. With ten stability deriva- 
tives all acting as variables in the inverse problem, a great deal of 
difficulty had been anticipated in effecting a convergence when using 
what is basically a simple iteration process, but to the degree of accur- 
acy obtained in the results mentioned above, no significant problems 
were realized. 

In future applications of this method, the procedures described in 
the sections "Development of Method of Analysis» and " Preliminary 
Test of Feasibility of Method’ would be useful to an investigator from 
two standpoints: they would provide data on the sensitivity of the specific 
configuration under investigation to changes in individual derivatives, 
and would provide valuable experience in manipulating the coefficients 


to obtain the desired response. 


Application of Method to Flight Test Data 


The four control deflections which were selected and the corres- 
ponding responses of the aircraft in the three variables : Ww , and 


D are shown in Figs. 15-18. As mentioned previously, a transparent 





overlay was made of each of these four sets of data to facilitate com- 
parison of the computer response with the flight data. 

One of the important features of the instrumentation was the fact 
that signals proportional to the four control deflections were available 
for forcing the computer by playing the recording of flight data, stored 
on the magnetic tape, through the decoding component of the telemeter 
ground station. From the output of the telemeter, the circuity shown 
in Fig. 13 was used to bias the steady state signals of the controls to 
zero and to provide proper scaling of the control deflection signal 
magnitudes. Without this feature, it would have been necessary to 
obtain the forcing function for the computer from an electronic function 
generator, with limited prospects for exact duplication of the aircraft 
control deflections. The use of the exact forcing function from the flight 
data thus eliminated a possible source of error by maximizing the com- 
patibility of the data collection system and the analog computer. 

The relationship between the actual aircraft response and the 
response of the analog containing the initial theoretical derivatives, 
when each was forced by a rudder pulse, is shown in Fig. 14. With this 
initial analog response (Fig. 14) as a starting point, the attempt to 
curve-fit the analog to the flight data was commenced. 

The initial steps of the curve fitting process were the same as the 
first four parts of the procedure previously described in "Development 


of Method of Analysis. Reference to Fig. 14 indicated that the most 
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obvious discrepancy was in the period of the Dutch roll oscillation. 
Since Tables IIl-a and I-b showed C no to be the derivative most 
effective in adjusting the period, the value of C n@ was increased 
by small increments until a relatively good period match was obtained. 
The values of the rudder control derivatives (C ns ; Cue , ye ) 
were then corrected to reflect the change in “ng » The next steps 
were to attack successively the discrepancies in damping, yaw-rate 
amplitude, and roll-rate amplitude by adjusting the appropriate deriva- 
tives in accordance with the information in Tables UlI-a and III-b. 
This procedure improved the match between the analog and the aircraft 
responses appreciably. However, in an attempt to further improve the 
matching, and at the same time gain further insight into the interplay 
of the effects of the various derivatives, all derivatives not adjusted 
in the foregoing procedures were successively altered to observe their 
effect on the analog responses. When no change, or an adverse change, 
in the responses was encountered, the derivative under investigation 
was returned to its original value. Those changes which appeared 
beneficial were retained. In each case, the effect of altering a deriva- 
tive, and the extent of the change were recorded. 

After a reasonably close match between the analog and aircraft 
responses to a rudder pulse was obtained, the analog was next forced 
by a rudder doublet, with the computer potentiometers still set on the 


derivatives established by matching responses to the rudder pulse. 
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The procedure which was then followed to minimize the discrepancies 
between the analog and aircraft responses was identical to that used 
for response to a rudder pulse, and again produced an improvement in 
the analog responses. 

The responses to an aileron doublet and aileron step function 
were then investigated in a like manner; that is, by first concentrating 
on the derivative changes indicated in Tables II-a and III-b, and then 
investigating the individual effects of all other derivatives. During the 
initial analysis of the aileron data, it was determined that an error 
existed in the aileron signal due to cable stretch in the control system. 
Corrections required for this error are discussed in Appendix C. 

During the conduct of the actual flight tests, it had been noted 
that in response to an aileron step function, the aircraft developed a 
nose-down attitude within two seconds after the application of the 
aileron displacement. This was caused by the large roll angle that 
developed. The resulting errors in the 2 and y traces were 
readily apparent and consequently, these two variables were disregarded 
in this case. The only matching that was attempted with aileron step 
forcing functions were the amplitude of the initial peak of the roll rate 
response, and the neutral stability of the spiral mode as indicated by 
the constant mean value of roll rate after the roll rate response had 


reached its initial peak value. The rapid development of the nose-down 


attitude when the aircraft was forced by an aileron step function could 
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have been avoided by the use of an initial condition, such as a steady 
sideslip. This would have extended the time during which the airplane 
was in the linear approximation range. Since the equations of motion 
of the analog are developed for perturbations from steady state, a 
similar initial condition would not be required when forcing the 
computer. 

After completion of one series of iterations for the values of the 
stability derivatives, the procedure was repeated for all forcing func- 
tions, again starting with the rudder pulse. It was noted that good 
matching was quickly obtained for any one of the forcing functions 
taken individually, but a shift to the other forcing functions offered a 
relatively poor match in the early stages. Gradually, however, as a 
result of the experience obtained, and by reference to the records of 
the effects of varying the individual derivatives, new procedures for 
establishing the value of some of the derivatives became apparent. 
These procedures represented a more definite method of evaluating 
the corresponding derivatives than did the influences indicated in 
Tables OI-a and UI-b, because in each case a particular response 
characteristic was controlled by an individual derivative. The deriva- 


tives and their individual effects are as follows: 


Ch - The period of the Dutch roll mode is controlled by 


8 


this derivative. This derivative has more than ten 


times the effect on the period than does any of the 


other derivatives. 






C, has a controlling effect on the damping of the 


r 


Dutch roll mode when the analog is forced by a 
rudder deflection. 
Changing this derivative changes the effectiveness of 


the rudder in producing yaw. Therefore, n 6, 

can be established by matching & and y ampli- 
tudes due to a rudder deflection. 

The magnitude of this derivative controls the match- 
ing of gb response in two respects; phase and ampli- 
tude. Phase error in the initial g response of the 
analog is affected considerably by this derivative. 
When the airplane is forced with an aileron step func- 
tion, the Dutch roll mode is excited to such a small 
degree that the initial peak of the D response can 
be attributed to the rolling mode and the magnitude 
of Clg. . ‘Fora given C1, : <i can be 
determined by matching the initial peak of the go 
response. Determination of the precise value of 
NS is explained below. 


Just as Cn. controls the damping of the Dutch roll 


mode, Ee controls the damping of the rolling mode. 


With Cl adjusted as described above, Cus can 


be determined by matching the o response to an 
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aileron doublet, especially during the latter half 
of the aileron input. 

The effect of this derivative is readily apparent in 
the gb response as an initial adverse gD due 
to a rudder deflection. Matching the initial form of 
the g response will thus establish le : 

In the response of the computer analog to a rudder 


deflection, ae has only a small effect. However, 


an upper and lower limit can be established for cae 


a 


beyond which the character of the Y trace is 
noticeably altered in magnitude. In the case of an 


aileron deflection, the effects of C and “nag 
a, 


Mp 


are so similar that one can be replaced by the other 

to a large extent. If an iteration is performed between 

these two derivatives, however, maintaining ae 

within the limits developed during a rudder deflection, 

until the best match is obtained in y due to an 

aileron doublet, C and C can be determined 
» Ynp ng, 

within reasonable limits of accuracy. 

Only a very rough value of this derivative can be 

established. Gross changes in C}.. will effect the 


stability of the spiral mode. The D response to 


an aileron step function is an indication of the spiral 
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mode stability: a constant mean value of , 
indicates neutral stability, a decreasing mean value 
indicates positive stability and an increasing mean 
value indicates negative stability of the spiral mode. 
Since the other three derivatives which affect the 
spiral mode (Cng j Cn,. Cle ) can be determined 
by the methods outlined above, the value of Cj. can 
be adjusted to match the character of the mean 
response to an aileron step function. 


C 
y 
é, 


"A 


- Because of the insensitivity of the responses to 


z changes in these derivatives, no refinement of the 


wee” ee eee 


theoretically determined values was possible. 


The above relations are summarized in Table IV. 


By the time sufficient information had been accumulated to 
develop the procedures outlined above, the initial discrepancies between 
the analog and the aircraft responses, for all four forcing functions, 
had been appreciably reduced. Subsequent use of the iteration proced- 
ure, in the case of each forcing function and its related responses, was 
confined primarily to adjustments in the derivatives which were based 
on the relationships of Table IV. In addition, when working with the 
responses to the rudder forcing functions, care was taken not to make 


large changes in those derivatives which had been established primarily 
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by matching responses to aileron forcing functions. This restriction 
was found necessary, since alteration of the aforementioned derivatives 
generally resulted in only minor improvements in the rudder cases, but 
had an appreciable adverse affect on the matching of the responses in 
the aileron cases. The same restriction was applied when matching 
responses to aileron forcing functions; inthis case, no extensive 
changes were made in those derivatives which had been established 
primarily by matching responses to rudder forcing functions. 

When a set of derivatives was found such that the responses of 
the computer gave satisfactory matching with the flight data for all 
four of the forcing functions, and no further improvement seemed justi- 
fiable in relation to the amount of time it would involve, the final value 
of each derivative was read off of the computer potentiometers and 


recorded. 
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RESULTS 


The final values of the lateral stability derivatives of the NAVION 
as determined by the method of this report are tabulated and compared 
with the initial theoretical values in Table V. The matching of the 
responses of the final computer configuration to the flight test data for 
all four types of forcing functions are presented in Figs. 15-18. 

During the iteration process, a better match had been obtained to 
the responses for each of the four forcing functions taken individually, 
but no single set of derivatives gave as good a match to the responses 
of all four forcing functions taken as a group as did the final set of 
derivatives listed in Table V. The major discrepancies remaining in 
the final matches are the generally low response of the analog in the 
Inns variable for all forcing functions; the too large amplitude re- 
sponse in the y variable for rudder doublet forcing functions; and 
the initial y response to the aileron doublet. It is felt that had 
more time been available, further investigation into the interrelated 
effects of the cross derivatives ~ and C,, would have resulted 
in nearly perfect response matches in the y and ra variables 
when using the rudder forcing functions. 

The accuracy of the final value of each derivative is directly 
related to the sensitivity of the system response to changes in the 


derivative. The derivatives C » GH Ge, Ce » Cay; 
Wis by A by A 
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Cj. , and C, are judged to be very accurately determined by the 
p i 
method. (See Appendix C concerning the accuracy of the specific 


values for Cl. and Di ayy of this report.) The values for Cl,» 


Cy , and GC, are susceptible to somewhat larger error due to 
"P OL, 
the smaller effects they produce on the response, and also due to 


overlapping effects in the case of Cnp and Cy Shame Finally, the 
a, 


value of the derivatives Cy and Cy._—- could not be established 

© ow 
with satisfactory accuracy. Other types of flight maneuvers might be 
advantageous in establishing these latter derivatives. 

The proposed method of extracting stability derivatives from 
transient flight data was determined to be completely feasible. As 
mentioned previously, there are inherent difficulties in the process of 
extracting all of the derivatives from flight data regardless of the 
method used. For certain applications, however, no better criterion 
could be established for obtaining the values of the stability deriva- 
tives than to match the response of the analog containing these deriva- 
tives to the response of the actual aircraft. In the field of servomech- 
anisms, for example, this criterion would be considered to be of 


prime importance. 


Another advantage of the method would be its economy in time, 
both in obtaining flight data and in extracting the derivatives from that 
data. In obtaining the flight data, any arbitrary forcing function can be 
used, providing it offers sufficient excitation of the modes of aircraft 


motion. Use of the telemeter data collection system permitted 
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convenient monitoring of flight data during the flights, and thus when 
necessary, instructions could be given to the pilot, to improve the 
responses obtained. Normally this would allow the recording of all 
the data for one flight condition in a matter of minutes. The ease of 
extracting the derivatives would tend to increase with practice. Un- 
doubtedly, Table IV, which contains procedures for establishing 
several of the derivatives, could be expanded with the aid of more 
time and experience with the method. For this reason, it is felt that 
better matching than is illustrated in this report would definitely be 
possible. It is important to note that the results of this investigation 
were obtained for the aircraft in the cruise condition only. Because of 
the large variation of certain stability derivatives with changes in 
flight conditions, further investigation of the method throughout a 
range of flight conditions is recommended. There is virtually no 
chance of computational errors affecting the results, for once the analog 
has been completed the computations are performed by the computer. 
Although the theoretical estimates of the derivatives were made 
as accurately as possible in this investigation, it was found that the 
rate of convergence of the iteration process was sufficiently rapid in 
the first stages to suggest that very rough estimates of the derivative 
values would be sufficient for use as initial computer settings. The 
time involved in refining the rough estimates would more than likely 
exceed the amount of computer time necessary to accomplish the same 


degree of refinement. 





Certain limitations of the method are readily apparent. The 
exact form of the equations of motion must be known in order to con- 
struct the analog. The application of the method to more complex 
aircraft would require further investigation and refinement; for ex- 
ample, any nonlinear effects of the aircraft system would have to be 
duplicated by electronic circuitry in the analog. The accuracy obtain- 
able for those derivatives to which the response is insensitive is 
undoubtedly of low order, but the possibility does exist of using differ- 
ent types of flight tests to establish these values. If these tests 
were compatible with the proposed data collection and reduction 


method, the additional effort could be minimized. 


CONCLUSIONS AND RECOMMENDATIONS 


The proposed method of extracting stability derivatives from 
transient flight data was proved to be feasible within the scope of this 
investigation. Matching of responses between the test aircraft and the 
equations of motion containing the derivatives is an excellent criterion 
for many applications of the derivatives. Due to the compatibility of 
the data collection system and the analog computer, the method is rela- 
tively fast and free from computational error. Only rough approxima- 
tions to the values of the stability derivatives are needed to construct 
the analog. In general the accuracy with which each of the derivatives 
is obtained is associated with the importance of that stability derivative 
in influencing the motion of the aircraft during the test. The use of 
four different types of control deflections as forcing functions is con- 
sidered a minimum requirement for accurate determination of all 
lateral derivatives. Subject to these restrictions, the exact form of 
each particular control deflection is arbitrary. Application of this 
method should include the preliminary procedures described under 


‘and " Preliminary Test of 


"Development of Method of Analysis ' 
Feasibility of Method. " 

Application of this method to non-linear aircraft would require 
further investigation and refinement. Any non-linear effects in the 


aircraft system would have to be duplicated by electronic circuitry in 


the analog. 


DA 


As a result of the investigation, the following recommendations 


are proposed: 


Ls 


Investigate the method at other flight conditions and with more 


complex aircraft. 


. Investigate the use of other control deflections to obtain 


responses which are more strongly influenced by the less 
sensitive derivatives. As an example, C1. might be more 
accurately obtained from an analysis of the response to a 


rudder step function. 


. Apply an initial condition such as a steady sideslip, prior to 


forcing the airplane with an aileron step function. 


. Augment the data (Table WI) which indicates the effect(s) of 


the individual derivatives on the response of the system by 
conducting a similar variation of the derivatives when forcing 
the computer problem with an aileron function. 

Improve the test aircraft aileron instrumentation by recording 
both left and right aileron angles and converting to average 


aileron angle by means of analog computer scaling circuits. 


yA 
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TABLE II 


PHYSICAL CHARACTERISTICS GF NAVION 


Ae WING 


1. Total Wing Area, Sy (includes flaps, ailerons 
and 19.87 Ft*+ covered by the fuselage) 


20 Span by 
3- MAC 


l. Angle of Incidence 
Root, iy 
Tip, dt 


5S. Twist 
Aerodynamic 
Geometric 

6. Airfoil Section 
Root 
Tip 

To Aspect Ratio, AR, 

9. Dibedral 

10. Root Chord 


lio Tip Chord 


A-1. AILERQN (one Aileron only) 
(Beaded skin, Frise Nose Balance, 


No Trailing Edge Extension) 
1. area, Sa 

2. Span, bg 

3. Aileron Deflection, 9, 
hh. Boost 

5. Trim Tab 


184.337 Ft.* 
33.378 Ft. 
68.352 in. 


+20 
~)0 


2° 311 
30 


NACA 15k 
NACA 6410R 
6 .0h4, 
0.5265 

705° 
Tod “Ft. 
3092 Fte 


120°-120° 


Wheel Throw 
Fixed bend tab on 
right aileron. 


2.161 Ft.% 
61.987 in. 
30° UP, 20° DN. 
None 


Fixed bend tab on 
right ailerm. 





Be 


Bol. 


Bo-2. 


C. 


HORIZONTAL TAIL 


l. Total Area 
(Includes 2.368 Ft.* covered by fuselage) 


2. Span, by 

3. Airfoil Section 

4. Taper Ratio, Az 

So MAC, Cy 

6. aspect Ratio, [eR ¢ 


HORIZONTAL STABILIZER 


le Area, Ss 


26 S y (with reference to FRP) 


ELEVATORS 


(No Trailing Edge Extensims ) 


lo Total Area, Se 

2. Spans De 

3. Deflection, Ge 

ho Trim Tabs (32" Span, h-1/2" Chard) 
5. Root Chard 


6. Tip Chord 


VERTICAL TAIL 


1. Total Area, Sy 


(Includes 2.577 Ft.* blanketed by 
fuselage and excluding 1.),83 Ft.¢ 
of dorsal fin) 


2. Airfoil Section 
Root 
Tip 


43.051 Ft.¢ 


13.172 Fto 

NACA 0012 
67 

3.34 Ft. 
4.02 


28.953 Ft. 
-},0 

Smooth Skin 

Flat sided 

No trim bungee 

Balance spring 

1.098 Ft. 

73-582 in. 

30° UP, 20° DN. 

30° UP, 30° DN. 
1.5 Ft. 


1.0 Ft. 


12.925 Ft.¢ 


NACA 0013.2 Mod. 
NACA 0012-6), * 


C-l. VERTICAL STABILIZER 
1. Area, Sy, 6.873 Ft.¢ 
(Includes 0.1427 Ft.* blanketed by fuselage) 
2.Ql (With reference to fuselage 7) 2° Nose left 
C-2. RUDDER Smooth Skin 
Rigged 3° Rt. to 
Fin g, Fixed 
Bend Tab. 
le Area, S, 6.052 Ft.¢ 
2. Rudder Deflection, & ; 17° L - 23°R 
30 Rudder Pedal Throw 5075 in. 
ho Trim Tab Fixed Bend Tab 
D. PROPELLER Hartzell 
lo. Activity Factor, A.F. 100 
2. Diameter, d 86" 
3e Pitch © @.75R 21° 30! 
Bo NOSE GEAR 
1. Travel + 30° 
2. Strut Angle (To Vertical) 10° 
F. MISCELLANEOUS 
lo Length, overall 27 025 Ft. 
2, Tail length, yy 16.88 Ft. 
G. POWER 
1. Engine Continental E-185 Engine 


185 HP @ 2300 RPM, 29" He 
M.P. @ Sea Level 





TABLE I0-a 


EFFECT OF VARIATION OF STABILITY DERIVATIVES 


ON AIRCRAFT RESPONSE") 


PARAMETER ee MAGNITUDE OF: “a 
INCREASE 


PERIOD 


slight 
decrease 





apne 


it oa, = en 


decrease decrease 












decrease strong 


decrease 











slight 
decrease 


ee ee a 

















strong 
decrease 


slight 
increase 


slight 
decrease 












moderate 
increase 


small 
decrease 


strong 
increase 


slight 
decrease 





















slight 
decrease 








decrease | strong 
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(1) Only rudder forcing function used. 





TABLE I[II-b 


STABILITY DERIVATIVES WHICH INFLUENCE 


SPECIFIC ASPECTS OF AIRCRAFT RESPONSE 


PEAK MAGNITUDE OF: PERIOD DAMPING 





Derivatives listed in descending order of effect. 


( ) -- Indicates a decrease in the variable affected. 


TABLE IV 


PRIMARY EFFECTS OF INDIVIDUAL STABILITY 


CONTROLLING 
DERIVATIVE 


“ne 


Cn, 


Cy 


e 
Plc 


“1, 


DERIVATIVES ON AIRCRAFT RESPONSE 


RESPONSE CHARACTERISTICS 


Period of the Dutch roll mode. 


Damping of the Dutch roll mode. 


Magnitude of WY and © response to a rudder 
deflection. 


Phase and amplitude of the $ response toa 
rudder deflection. 


Amplitude of initial $ response to an aileron 
step input. 


Having established a value for C1. > Clp 
can be determined by the roll respmse to an 
aileron doublet, particularly during the last half 
of the control input. 


Adverse roll characteristics, which are apparent 
inthe ¢ response immediately after a rudder 
deflection. 


An upper and lower limit for Cnp can be 
established by an amplitude match of the 
responses to a rudder deflection. Cyn, can then 
be adjusted within these limits along with the 
adjustments in Ch aw to provide the best match of 


y response to an aileron doublet. 


(continued) 





TABLE IV (continued) 


Si Response insensitive to this derivative, but a 
qualitive match of the slope of the mean 
response to an aileron step function is the best 
means available for determining Cj, . 


ye Response too insensitive to these derivatives 
to permit the determination of a special effect 
a» * in either case. 
y 





TABLE V 


THEORETICAL AND FINAL EXPERIMENTAL 


LATERAL STABILITY DERIVATIVES 






Final 
Experimental 






Theoretical 





Derivative 





* From Ref. 9 


= Dee aia ae, 





< 
‘eee | 
. SZ 


33° 4.59" 


‘ 4 | 
7) | 
PROP DIA 7 425 Ye he 


a 


p— 
t Ss it 


— . or + 
<— [OA GGA - 


FIG 1 
GENERAL THREE-VIEW DRAWING 
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APPENDIX A 


THEORETICAL AND EXPERIMENTAL DETERMINATION 
OF THE DYNAMIC RESPONSE CHARACTERISTICS 


OF THE SIDESLIP VANE 


General 

The instrument system chosen for measurement of the sideslip 
angle consisted of a light vane mounted on a boom which extended forward 
from the left wingtip. (See Fig. 4). As has been previously stated, it 
was considered of prime importance that the sideslip measurement sys- 
tem have dynamic response characteristics such that no coupling error 
would be introduced. 

Since no data were available for the response characteristics of 
the vane selected, theoretical calculations were made, and experimental 
tests conducted in order to determine if the natural frequency and damp- 


ing ratio of the vane were within acceptable limits. 


Theoretical Analysis 
(See Ref. 1, pp 11.51-11.55). Notation as per sketch below. 
Vv = true airspeed, ft/ sec 
Vj = indicated airspeed, ft/sec 
M m moment about pivot 
Ch = normal force coefficient 


iF = moment of inertia of vane 





Mon. 
| <<, 
lon. O on 
Veer ~SL06 jy Qe — Vane Pivot Axis 
a rP lane uy Disturbing Moment 
4 Xi § _—- Acrodynamic 
Normal Force, N 
Earth Reference ¥ 





AX!S ae Vane Lifting 
Surface. 
P SESS 


C 
Aerodynamic 


Chordwise Force 


(1) Oo = = a (~ However, due to angular velocity of the 


vane with respect to the air mass, the angle of attack at the vane, 


as os si r 


2) oe + 


V 
7 oe ae 
(4) Aes a C. N . , dynamic pressure 


es i 
\ aS } 4 G 
(5) a aa Ce N S 9 Sy, uty 
dX, 


6) T\ -4AN —-( 64 rae 


combining (3), (4), (5), & (6), 


vs f 








mol 4G 4 
a 


ae: \ 
sak he Cn oe _ Ge 
om 





where: ee mene een artnet 
& Se aa 
o ann =“ % — ao. } Natural Frequency 


of Vane 











— on aries Damping Ratio 
CF ee > / of Vane 
aime ae 2 ec. 
ado) \\ = yO SS iG \] L ZX — Aerodynamic 
aes Stiffness of Vane 
Conditions for Analysis: 
Density Altitude: 5000 ft 
V;: 120 mph = 176 ft/sec 
y 189.6 ft/sec 
S 0361 £t* 
He 157 ft 
AC J . 
sae 1.72/radian (See Ref. 5) 
I : 6.33 x 107° slug £2 (Determined experimentally) 
From (8) above, 
proenecncentntecncnernmmemn tennant ne 
, oC ,,.002578 .* 036) —oisis =~ 1672 
2 | 2~« 6.33 x 10-5 
\ 
cn). = 5 oe 
sec 
LO 
ee = é =e 12.01 Cycles/se¢ 
vx a T cae ee 
From (9) above, 


aS | .002378 4 , 8616 y , 036] « 157)? 2 1.72 
a 8 . 6.33 


4 107 
\ s 031 














Experimental Analysis 
Since the above theoretical dynamic analysis of the response char - 


acteristics of the sideslip vane was based on a moment of inertia of 
limited accuracy, and neglected the effects of friction damping in the 
system, it was considered desirable to determine the dynamic response 
of the vane experimentally. Ref. 1 (p 11.55) reported that wind tunnel 
calibration of vanes is usually unsatisfactory because of the poorly- 
damped vane response to disturbances, and the usual high level of tur- 
bulence in wind tunnels. For this reason, the sideslip vane system was 
modified to permit in-flight dynamic testing of the vane when installed 
on the airplane. <A spring and plunger mechanism, actuated by a sole- 
noid, was installed on the boom, as indicated in Fig. 4. When in flight, 
operating at the desired test configuration, the vane was deflected by 
sideslipping the aircraft, and locked from the cockpit by electrically 
extending the plunger. (See Solenoid Switch, Fig. 3.) Thus, when the 
solenoid was reenergized, the vane was released, and responded as if 
to a step function. Fig. A-l shows the transient response of the vane to 
such a step input. Calculations for the experimental determination of 
the natural frequency and damping ratio of the vane from the transient 


response of Fig. A-l areas follows: 


Conditions for Test 


Density Altitude: 5000 ft 


Vi: 120 mph = 176 ft/sec 





In Ref. 6, it is shown that the undamped natural frequency and the damp- 


ing ratio may be determined from transient data as follows: 


a | 
a aa 


D 











t= 
e. 
a c= i: 
W = = aN 
5 = 
ae 
where 
> yk ioe. “epee ais t 
XL pari Cc ecTremen 


From Fig. A-l, 


| _ .065 8ec., Xy = 4. 0, Xo = 3.1 


4.0 
~ = BCA | aad 2546 





| ae 2 
(Sy = 2st + TT. 97.0 Ra/cec. 


\ .065/2 














{ - é = 15.44 Cycles/ 
ry ane sec, 
O a 2546 
ee = i a .08 
>) _ 2h, ~~} 3.15 _- 081 





Comparison of Results 


| Theoretical | Flight | 

Analysis | Test : 

fn, Natural | | : 
Frequency | 12.01 | 15.44 

cps | | : 

| | 

~—<e | 
Skeet 0.031 | 0.081 | 
Ratio ! | | 

Ee 


t 
P] 
aE A I OS EO OS BY GF hg TI et fh ere hy TE IR my eS ee NE Ee 


} 
{ 
) 
! 
‘ 


Concluding Remarks 


(1) The natural frequency of the sideslip vane has been shown to 
be of the order of 16 cycles per second, and the damping ratio 
approximately 0.08. | 

(2) The ratio of the natural frequencies of the sideslip vane, and 
the airplane in the lateral oscillation is approximately 40. 

(3) In view of (1) and (2) above, no dynamic correction to sideslip 


vane data was required. 
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APPENDIX B 


DETERMINATION OF LATERAL STABILITY DERIVATIVES 


FOR THE NORTH AMERICAN NAVION FROM THEORETICAL 
FORMULAE, AND FROM WIND TUNNEL AND FLIGHT TESTS 


FROM GEOMETRICALLY SIMILAR AIRPLANES 


The determination of stability derivatives for an aircraft from 
theoretical and experimental data is a problem which has received wide- 
spread attention. Numerous publications, outlining a considerable num- 
ber of methods for calculating or estimating derivatives, and presenting 
the results of extensive test programs, have appeared over the past 
twenty years. (See bibliography to Ref. 7). However, in spite of such 
a wealth of information on methods and test data, estimation methods 
are often found in practice to yield only fairly accurate values, suitable 
for making only first approximations to dynamic stability characteris - 
tics. In addition, the application of several theoretical methods to the 
Same airplane often results in wide discrepancies between some of the 
calculated derivatives. 

In general, the most reliable methods for calculation of lateral 
stability derivatives depend on wind tunnel data for accurate values of 


the sideslip derivatives, Cy © . Cra es) C . This is par- 


) 
ticularly true for propeller driven aircraft, where the propeller lateral 


force and slipstream effects are important contributions to the sideslip 


derivatives. No satisfactory purely theoretical method has yet been 





developed for accurately estimating these derivatives for a complete air- 
plane primarily because of large interference effects between the various 
airplane components. 

The remarks above concerning wide discrepancies were found to 
be particularly pertinent to the calculation of the lateral stability deriva- 
tives of the North American NAVION. The derivatives for this airplane, 
including the sideslip derivatives, were calculated by five well known 
methods, and the results compared with test data, both for the NAVION, 
and also for several geometrically similar airplanes. The results, pre- 
sented in Table B-IL, showed a variation of at least 100% between the high 
and low values in the case of all derivatives with the exception of CG P 
which varied by approximately 30%. 

It is to be noted in Table B-I that the method of Ref. 5 gave the 
best correspondence between theoretically calculated derivatives, and 
the "average" derivatives of the several methods and tests results. 
Therefore, the detailed calculations for the lateral stability derivatives 
of the NAVION by the method of Ref. 5 are presented in the following 


pages. 





GENERAL PARAMETERS FOR NORTH AMERICAN 


NAVION N91566 


In addition to the physical characteristics of the NAVION pre- 
sented in Table II to the basic report, the following parameters, 
determined from the assumed test conditions, are required for cal- 


culation of stability derivatives. 


W ..... Gross weight during test, lbe--+--ececese 2770 


hg----- Density altitude, ft ........ SACO 0: ee 5000 


e .... Air density, slug/ ft? ...4) eee ee «002049 


Vow acncated airspeed, mph ..,ce wea. -.: eee TAt 
V ..... True airspeed, ft/sec o-sscecsceceserece 189.6 
Cr .... 2 otal lift coefficient ....... ee ee ee 0.407 


Nr . .. Tail efficiency factor (estimated) ....... 0.90 


Unless otherwise specified, all further references in this appendix 


are to Ref. 5 to the basic report. 


TABLE B-II (a) 


THEQRETICAL DETERMINATION OF THE LATERAL STABILITY 
DERIVATIVES GF NAVION N91566 


“= 
| (5) (8) 
: Reference Perkins § USNPOGS 


| & Hage | Notes 


; 









i 
| 
| 
j 
{ 
a re. 
{ 
‘ ' | 
{ 
| 





| | 
| | -.592 = -.612 
} 
061 a <0 
: oe CG : 3 | 0708 - 0409 a 0057 
| 
| Ce e | ~.0816 | ~21236 | -.0981 | -,1137 : -00991 | =.107 
CO | ‘ : | 
| Re | eG 495 | = ~.6 ~.30) am 
{ j 
| Cr | e102 ollO = e102 | 00978 | e103 
| 
| — : -,051 : -.02);3 is -.0431 | ~,022h -.035 
| i : 
ena. | ~.0937 : -.0945 = -.0953 | = -.095 
Pog | | : | 
Pe ° = = | 
| DA | ote ! a i e111 
Cn iL | : 
Qu. | 00087 | = oa aed H = - .0087 
| 
Bos : 
| -.0 = = 
Sn, | 790561 - : - : -.0561 
C 133 - - 
Me | - /@ 


* From Ref. 8. 








B=5 


TABLE B II (b) 


COMPARISON OF LATERAL STABILITY DERIVATIVES 
FQR AIRCRAFT GEOMETRICALLY SIMILAR 
TO NAVION N91566 





ee -—— 
) 





a ee an ere 8 er 0 te 0m ee ce pee we 


(13) (4) (15) 
™ RM 












16 
REFERENCE P.Re* TN i 
232 2195 1327 L53K09 
Co | 4390 | 633k 40 133 
>"/Su | 068 | .088 | was waza | 0 
Ave/ | 
506 e460 ol hyl 156 51 
ae | 60h | 6.25 | 5.9 | 5.89 6.25 
x | 527 
| [eee | 7.5 | 
0 
ey, fo 





STABILITY DERIVATIVES 


——— Se a nS oe ae AR at Am I SE RE del fhe 


| 
: 
: 





Oye | -.595 | -.802 -.573 | -.515 

or e 141 »132 20975 01146 

CR -.0822 | -.08 | -.0688 | -.0573 

: Re -243 | | -39 

Se 203 i 06 

Cre | ~2007 -.02 | -.020 -.018 
Crp, -.160 13 | -.138 ~ 2143 
Opes | «203 | 0203 
Cale : .0166 0166 
Oke : ~.107 -.107 
Cm | 101 : | 101 
Cr - : eel! : | =.0035 


* For Navion 5113 





oye 
NOTE: Ref. 5 does not present a theoretical method for calculation of 


Cx C . Therefore, sare was calculated from the method of 
Ref. 8. 


fee — 126 L 
re i S . where 


| «= Fuselage length, Ft, . ...... . 28.25 


GY _ zk2 (33.38) (27.25) | _. 592 


& 184.2 


le 


= 


Total directional stability, “ne » is composed of contributions 
from the wing, fuselage, propeller, vertical tail, and interference cor- 
rections due to wing configuration and to sidewash from the fuselage. 
In equation form, 


(Cre) 


airplane 


CC. c) anes Care + Cao.) Cre) es A G ve + A Lae 


prop* 


(ay Cy . -.0006 | A| 2. 57.3 =0 
‘ew 
(b) A Cn = 0 for low wing aircraft (see p. 320) 


(c) CA fuselage: An empirical relationship for this contribu- 
tion, which is particularly applicable to single engined aircraft, 


is as follows: 
2 (we? 
S L W 
=» .96K (22\( £ (2 (may 
nee fus. eC 3 abr L2 Wi 


where 


Ss = projected side area of fuselage in Ft? 


Le = fuselage length 


ly, L2; Wh, W2 from accompanying sketch 


Ke from Fig. 8-4, p. 320 








For the Navion, the above quantities are as follows: 
i _)  Sra 2, - - 
Lig, ee rrr |S 

Ly, en ee ee ee ae 4.25 


ft 
ft 

OE ee ee eee es 
es en >. oOo 
ft 


Cr. + + slew cove ge Le. 1.76 


d5, nose toc.g., ft. ..,.,..,4.4. 8.05 


Ko, (from Fig. 8-4). --......... o 10 
therefore, ; ; 
L 
C = -(,96) (.110)/ 96.22 (27.25 (4-28) ( iS) 
“@ fuse nee ) (3528 33. 38/\ 2.61 ZIT 
eG . 7+ 0434 


"G fuse 
(d) (no prop from pp. 321, 


" ACY 
a ° i¢ 
aoe fe : rai ay P ¥ > 


5st 


D = Propeller diameter, ft oe © ee © @ @ © «@ 7.0 
N = Number of propellers ......... 12 


4 


= Distance from c.g. to propeller plane, ft, 7.90 


~) 


‘- 





(ec c) = 4. ( =) 
de ~ aoe ite where, from p. 321, 


© 


k «1.5 for full power, 


cA — g 
“a « .00165 for two bladed propeller 


For Navion N91566 at full power, 


ene = 185, “].= .632 (Ref. 9) 





~T a 550” SH 550 x .632 x 185 
——"7 . 00204 189.6 3 x 7 
CIE 1s ae © 


G full power 0940 


For cruise power, 


SHY = 130, “We = 18 
550 x .78 x 130 
ay ~ 5002049 x (189.6)3 x (.7)% — adcls 


K cruise power . le cruise power K (1.5) 


a full power 


- 0815 
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(e) Cre vertical tail 
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OCu..4, the slope of the side force curve for the vertical tail, is a 


function of the effective aspect ratio of the vertical tail, Be e! 


Z 
RQ e * 1-55 byt where Sy,t is the nominal vertical tailarea 
Svt 


a3 indicated on p. 325. 
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From Fig. 8-8, p. 324, 
Covt = -048/deg = 2.75/radian 
X yt = distance from airplane c.g. to aerodynamic center of 
the vertical tail. 
As measured from a scale drawing for appropriate c.g., 


= 16.88 Ft. 
X vtNavion 


“ ve * ee , assumed = .90 
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When sidewash effects are included as a part of the tail efficiency 
factor instead of as an increment, A 1 Gr C , the total 


efficiency factor can exceed unity, and is of the form: 


8 See for example, Fig. 2la of Ref. 10. 
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(f) A 2 re » sidewash or interference flow correction due to 
wing-fuselage combination, 


From Fig. 8-7, p. 325, with Dg and d, from part (Cc) above: 
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where 


Kae pn. is a correction for wing-fuselage interference due to wing 
\ \~ 
height. 


A & “te is a correction for wing-vertical tail interference as a 


function of wing height. 
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The effects of power on “ j , are negligible for the "clean" 
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condition. 
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~ Yeawing ‘ ao} = -.0968/ra, 


N C Xe tipshape “ % from Fig. 9-5, p. 345. 


(b) \ Ske : From Eqn 9-6, p. 347, 
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po 4 C low wing 
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Cs (~ vertical tail 
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where 2,4 «= distance from X axis to aerodynamic center 


of vertical tail. 
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(d) INS ‘a : From Eqn. 9-H, p. 346: 
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Co Gs » Damping in Roll 


See Fig. 9-14, p. 351. 
For A =0.5, Ae = 6,04, 
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By strip integration, and assuming elliptical lift distribution, 
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elliptical lift distribution. 
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M7, Damping in Yaw 
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where AS = rudder effectiveness, 
Weom Fig. 5-33, p. 250 
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oe = aileron effectiveness, a function of aileron-wing 


chord ratios. 


(S4 « \ is a function of the location of the inboard and 
“No. outboard ends of the aileron in percent of 
semispan. (See Fig. 9-3, p. 344.) 
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From Ref. 3, for the Navion, 
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APPENDIX C 


CORRECTIONS TO RECORDED AILERON DATA 


During preliminary matching of responses of the test aircraft 
and the computer to aileron forcing functions, it was determined that 
the recorded values of aileron position angle would require correction 
for errors introduced as a result of a combination of aileron cable 
stretch and a non-linear hinge moment curve. In order to quantita- 
tively establish the extent of the required correction, the computer 
problem was forced by an aileron step input to the left and the roll 
rate response matched to the flight data by adjusting a 
a Similar aileron input to the right was then introduced, and was 
observed to cause a computer response in roll rate that was smaller 
than the flight data by a factor of 0.78. This factor of 0.78 was 
proved to remain constant over a wide range of values of Clpy ; 
the one stability derivative that would have a marked effect on the 
roll rate due to aileron deflection. Similarly, the roll response to 
an aileron doublet could be brought into perfect match with the flight 
data during the initial roll to the left, but no reasonable value of 
oT would produce as high a roll rate to the right as was recorded 
in the flight data. (Values of ois as low as 0.28 were investigated. ) 


At this point, equation (9-33) of Ref. 5 was rewritten in ratio 


form, and applied to responses of actual flight data of aileron step 





functions to the right and left, as follows: 


D tet y. da, left 
D right 5a right 


The k factor was inserted in the equation knowing that its value 
should be unity. However, upon solving for k it was again deter- 
mined to be the factor 0.78. 

This analysis offered sufficient proof that an instrumentation 
error had been included in the data. The gyro measuring the roll 
rate response was known to be highly accurate, but aileron deflec- 
tion was measured directly at the right aileron and the calibration 
curve was obtained on the ground under a no-load condition. Aileron 
cable stretch was known to exist to some degree in the aircraft used, 
and for this reason, aileron deflection was measured directly at the 
right aileron. Instead of alleviating the instrumentation error, how- 
ever, this merely caused it to be additive during left roll and sub- 
tractive during right roll. This condition can be demonstrated by 
first considering the non-linear hinge moment characteristic of a 


modified frise aileron: 





It should be noted that the change in hinge moment is nearly pro- 
portional to the deflection when the aileron is displaced down. There- 
fore, due to cable stretch, the down-going aileron would be attenuated 
by a constant factor in comparison to the no-load calibration. The up- 
going aileron has a relatively constant hinge moment over the range 
used in this investigation and would therefore parallel the no-load cali- 
bration curve. The resulting error is illustrated by a plot of control 


deflection vs. aileron deflection: 
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When the ailerons are deflected to the point where the right aileron is 
at position A, the signal output is proportional to the length 1-3, but 
the actual deflection is the shorter length 1-2. Correspondingly, if 
the controls are reversed to the point where the right aileron is at 

point B, the signal output is proportional to the length 6-7, but the 


true position is then the larger value 4-5. 


Qualitatively, this discussion would hold true whenever the slope 
of the hinge moment curve differs between the up-going aileron and the 
down-going aileron. Although no hinge-moment data was available for 
the Navion, the ailerons are of the modified frise type, which are 
known to have this characteristic. 

Having qualitatively explained the error in terms of non-linear 
hinge moments, and having quantitatively found the magnitude of the 
error both from flight data analysis and from the computer responses 
to aileron step functions, it then became necessary to remove the 
error from the signal in order to continue the computer analysis. 
Although the ailerons had a differential gearing characteristic for large 
deflections which would affect the apportionment of the error between 
the right and left Signals, the differential effect was less than ten per - 
cent in the small range of deflections used, and thus the apportionment 
of the error could be made on a 50-50 basis with reasonable assurance. 
(This concept could be verified by making reference to the plot of con- 
trol deflection vs. aileron deflection. With no differential in the 
aileron system, the percent of error for a given control deflection to 
the left is equal to the percent of error for a corresponding control 
deflection to the right, but the signs of the errors are Opposite. Adding 
a differential pearing effect changes this equality between the error 
signals, but not to a marked degree.) 


To effect the change in the aileron Signal, circuitry was provided 


to divide the signal into two paths, one path for the signal during right 
roll, and another path for the signal during left roll. The right roll 
signal was then attenuated by the previously established factor of 0.78. 
This circuitry is shown in Fig. 9. 

With this change the character of the signal had been corrected, 
but the magnitude of the entire signal was known to be low because the 
correction had all been placed in one direction as an attenuation rather 
than being apportioned evenly between the two directions of roll as an 
attenuation of one and as amplification of the other. The amount by 
which the signal was low was the ratio of 0.78/0.89, or (0.877), for 
dividing the signal by that factor would halve the previous attenuation 
in the right signal and amplify the left signal correspondingly. Thus 


the correction would be apportioned equally to the two signals: 
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Instead of applying this second correction factor to the signal which 
would have required additional changes in the circuitry, the investiga- 
tion was carried to its conclusion and the correction factor was then 


inversely applied to the value of the aileron derivatives, the net effect 





of the two alternative methods of correction being identical in the 
final result. 

In conclusion, it should be stated that the corrections made to 
the signal were justified by the viewpoints brought out in the above 
discussion. However, to state that the error was completely elimin- 
ated would seem overly optimistic, especially with respect to the 
apportionment between attenuation of the signal in one direction and 
amplification of the signal in the other direction. It is strongly felt, 
however, that any error remaining in the aileron signal had no effect 
on the value of any of the derivatives other than the aileron derivatives, 
Cis and C, . The magnitude of the error produced in these 

a Oa, 
two derivatives should not exceed five to eight percent. 

A more suitable method of instrumentation would have been to 
include a transducer at each aileron, and sum the two aileron angle 
Signals and compute the average angle as a function of time on the 


analog computer. 




















